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We study the vortex state of a layered superconductor with vertical line nodes on its Fermi 
surface when a magnetic field is applied in the ab-plane direction. We rotate the magnetic field 
within the plane, and analyze the change of low-energy excitation spectrum. Our analysis is based 
on the microscopic Bogoliubov-de Gennes equation and a convenient approximate analytical method 
developped by Pesch and Dahm. Both methods give a consistent result. Near the upper critical field 
H C 2, we observe a larger zero-energy density of states(ZEDOS) when the magnetic field is applied 
in the nodal direction, while much below H C 2, larger ZEDOS is observed under a field in the anti- 
nodal direction. We give a natural interpretation to this crossover behavior in terms of momentum 
distribution of low-energy quasiparticles. We examine the recent field angle variation experiments 
of thermal conductivity and specific heat. Comparison with our results suggest that special care 
should be taken to derive the position of line nodes from the experimental data. Combining the 
experimental data of the specific heat and our analyses, we conclude that Sr2Ru04 has vertical line 
nodes in the direction of the a-axis and the b-axis. 
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I. INTRODUCTION 

Unconventional superconductors are one of the most 
important materials in the modern condensed matter 
physics, as a key to understanding strong electron cor- 
relation effect. A number of superconductors, high- 
T c cuprates, heavy-fermion metals, ruthenates, and or- 
ganic compounds exhibit unconventionality in the sense 
that the superconducting gap vanishes somewhere on the 
Fermi surface, resulting in power-law behaviors in vari- 
ous thermodynamic quantities. However, while many su- 
perconductors were found unconventional, detailed gap 
structures are still unexplored except for the cuprates. 
One of the difficulties in clarifying the gap structures 
seem to lie in the lack of experimental probes sensitive 
to quasiparticle momentum distribution. 

Recently the vortex states have been attracting much 
interest, because the positions of gap nodes can be de- 
tected. When a magnetic field is applied parallel to the 
superconducting plane(ab-plane), various physical quan- 
tities depend on the angle between the magnetic field 
and the superconducting gap nodes. Hence, by rotat- 
ing the field within the plane and tracing the change of 
physical quantities, one can obtain the information of 
the gap nodes. So far, thermal conductivitjiSiSiiiSiSiLSiS 
and specific healii * 1 ^* 1 ^ have been measured in the 
vortex state for a number of layered unconventional 
superconductors including Sr 2 RuQ 4 4 i 5 i 10 , CeCoIn j 6 ! 11 , 
K-(ET) 2 Cu(NCS)2i, YNi 2 B 2 C&ia and PrOs 4 Sbi2£. 

However, it is found that some of these experiments 
show the behaviors incompatible with the results of 



theoretical analyses. For example, microscopic cal- 
culations show an existence of vertical line nodes in 
Sr 2 RuO^*i^ii^*i&, while the experiments suggest line 
nodes run horizontally 4,5 . Probably these discrepancies 
are attributed to the lack of firm theoretical basis in an- 
alyzing the experimental data. Actually, experimental 
data have been interpreted based on a phenomenologi- 
cal Doppler-shift method, which has been claimed to be 
quite unreliable in some cases 1 ^. So far, there have been 
few microscopic analyses on quasiparticle states under a 
magnetic field parallel to the ab-plane. In particular, no 
microscopic analysis has been done on layered unconven- 
tional superconductors. Therefore, it is crucial to estab- 
lish a reliable theory in these systems, and give a correct 
interpretation of the experiments. 

In this paper, we present a detailed study of quasipar- 
ticle density of states in a layered superconductor under a 
magnetic field. We will focus on Sr 2 Ru04, in which posi- 
tions of the line nodes are still controversial. A cylindrical 
Fermi surface with vertical line nodes of f-wave symmetry 
is assumed, and a magnetic field is applied parallel to the 
ab-plane. We concentrate on the two cases where a mag- 
netic field is in the nodal direction, and in the antinodal 
direction. We investigate the low-energy quasiparticle 
states, and apply the results to the interpretation of the 
experimental data. 

Our analysis is based on the microscopic Bogoliubov de 
Gennes equation and an approximate analytical method 
invented by Pesch and recently developped by Dahm. 

In the next section, calculational formulation is de- 
scribed. In section ITTT1 we will show our results and dis- 
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cussion on the experiments. Section llVI contains conclu- 
sion. 



II. FORMULATION 

A. Model and Some features 

First, we introduce some general features of our model. 
We study a quasi-two-dimensional layered superconduc- 
tor which has a cylindrical Fermi surface with small c-axis 
dispersion. Thus we assume a dispersion relation, 



Pab 

2m n i 



v r cos k z 



(1) 



where p a b is the momentum in the ab-plane and the c- 
axis wave number k z varies in the interval [— 7T, 7t]. In this 
system, the Fermi velocity can be written in the following 
form, 



wp = vp(cos 4> e a + sin <f> + e sin k z e c ), 



(2) 



with Dp = and v c — evp. Here, the azimuthal angle 
<j> varies between [0,27r]. The c-axis dispersion e exists 
due to a small inter-layer hopping. 

In this paper, we consider a spin-triplet superconduct- 
ing order parameter with its d- vector directed parallel to 
the c-axis 



A(k) 



A(k) 
A(k) 



(3) 



We fix the momentum part of the order parameter as 

A(k) = A (k a + ih)(k 2 a-k 2 b ), (4) 

and ignore field-induced symmetry change of the order 
parameter. This is the simplest model of chiral state with 
four-fold symmetric vertical line nodes. In this model, 
nodes exist at \k a \ = \kb\- However, exact positions 
of nodes are not important in the following discussion. 
Qualitative behavior of low-energy density of states is 
determined by the relation between applied field and line 
nodes. 

In order to study the vortex state under a magnetic 
field parallel to the ab-plane, we assume a spatial vari- 
ation of the order parameter, ■0(r)A(k). Here ?/>(r) is 
described by the Abrikosov vortex square lattice with 
anisotropic superconducting coherence lengths, 
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where, £ c and £ a {, are the superconducting coherence 
lengths in the c-axis and the ab-plane directions, respec- 
tively, and v± denotes the coordinate for an axis which 
is in the ab-plane and is perpendicular to the magnetic 
field. For example, ip( r ) is equal to zero at the position 
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^£ c , representing the center of a vor- 
tex. In eq. (5), the prefactor 2* is a normalization factor 
to let the cell average of IVK 1 ")! 2 equal to 1. We have 
adopted the Landau gauge for the vector potential as 



A(r) = Br ± e c 



(6) 



It is apparent that ip(r) has a periodicity with respect 
to r, which corresponds to the vortex unit cell. The spa- 
tial periods L c (c-axis direction) and L a b (ab-plane direc- 
tion) are related to q as 



<?£afc, 



2tt£ c 



(7) 



(8) 



We choose q so that the ratio is equal to ^ i n t ne 
whole range of magnetic field. Then, tp(r) can be rewrit- 
ten as 



i/j(r) = 2 3 y. exp(2i7rn 7r(- — 

L * V L,. L„h 
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Since the size of a vortex unit cell is inversely propor- 
tional to the average induction B, we have the following 
relation 



Lj(B) = V2n^ 



B C 2 
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(j = c,ab). 



(10) 



Finally, we introduce parameters and physical quanti- 
ties of interest. Most properties in our system are deter- 
mined by a single parameter, namely the reduced order 
parameter rj{B) defined by 



2 A (B)L ab (B) 

IT Vp 



(11) 



where A (B) is a spatial average of the order parame- 
ter. Note that rj monotonically decreases as increasing 
magnetic field H and becomes when H reaches H C 2. 

We are interested in how the density of states(DOS) 
depend on the angle a between the field and the line 
node when a field rotates within the plane. It is expected 
that the density of states oscillates with the period ^, 
reflecting the fourfold symmetry of the order parameter. 
Therefore, we can concentrate on the two cases: when 
a field is applied in the nodal direction(a = 0), and in 
the antinodal direction (a = 7r/4). DOS is considered to 
take its minimum and maximum in one and the other of 
these two cases, respectively. 



B. Bogoliubov-de Gennes equation 

Using above dispersion relations and spatially inhomo- 
geneous superconducting order parameter, we solve the 
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Bogoliubov-de Gennes(BdG) equation, which is consid- 
ered to be the most reliable approach^. 



H u a (r) + A(-W) il>{ 



A*(iV) ip*( 



= Eu a (r), (12) 

, -H$ V<r (r) 

= Ev a (r), (13) 



where 



2m a b \ c J 



(14) 



Here we assumed e«l, and neglected the effect of the 
c-axis dispersion. This prescription corresponds to ne- 
glecting a coherence along the c-axis direction, or in the 
quasiclassical sense, to taking account of only the tra- 
jectories parallel to the ab-plane. We will discuss the 
details of this prescription later in the next section. Eq. 
(12) and (13) can be decoupled into 2x2 matrix equation 
for and (u^v^). Since both the pairs satisfy the 

same equation, we can work on only one of them, say, 

We numerically diagonalize Eqs. (11) and (12) by de- 
scretizing the coordinate r, and obtain sets of eigenval- 
ues Ek and eigenfunctions (uk{y),vk (r)). Using the ob- 
tained eigenfunctions, we calculate DOS vbcig{^)'- 

v BdG (e) = f dr[\u K (r)\ 2 S(e-E K ) 

rri ^ n J 

-\v K (r)\ 2 6(e + E K )]. (15) 



Ek>0' 



Due to the translational symmetry along the magnetic 
field, momentum parallel to the vortices, pu, becomes a 
good quantum number. Hence, for each cigenfunction, 
we can define an angle between the magnetic field and 
quasiparticle momentum as 



arctan ( 



P 2 F-Pl 
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(0 < 9 < tt), 



(16) 



where we limit the range of 9 to [0,7r] due to the reflec- 
tional symmetry about the magnetic field. Since momen- 
tum normal to the magnetic field p± is not a conserved 
quantity, p± has a finite width 6p± for each eigenstatc. 
Nevertheless, is much smaller than 1 (of the order 



of 



feirfa 



not too much below T c . Therefore, we can 



consider 9 as a well-defined quantity in the quasiclassical 
meaning. Using this 9, we discuss the momentum distri- 
bution of the quasiparticles contributing to zero-energy 
density of states. 



C. Approximation due to Pesch and Dahm 

Before showing our results of BdG equations, let us in- 
troduce an approximate analytical method, invented by 
Peschii, and developped by Dahmi& ? . We will compare 
the obtained results with this approximation. Near H C 2, 
spatial variation of order parameter is small. Hence, in 
the Eilenberger equations, it is allowed to replace normal 
component of quasiclassical Green functions g by its spa- 
tial average over a vortex unit cell. With this averaged 
quasiclassical Green function, we can calculate various 
observable quantities (e.g. DOS) in the averaged form 
over a vortex unit cell. According to Dahroi&, even much 
below H C 2, this approximation gives the result quantita- 
tively in agreement with that from the rigorous Eilen- 
berger equation. 

Here, we summarize the main results of this method. 
For details, see RefE We assume that a spatial 
variation of the order parameter is described by ip(r), 
Abrikosov vortex lattice introduced in the section HTK. 
Then, the averaged density of states vp£>{e){ya the unit 
of DOS for the normal state vq) can be written in the 
following form, 



fpc(e) = (Re 



1 



v/l + P(v F ,-. 
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where (• • ■ )f means averaging on the Fermi surface, and 
in our system, -P(v F , iuj n ) is written as 



4|A(k)| s 
7t|v f _l| 2 
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where vpi is the projection of the scaled Fermi velocity 
onto the plane normal to the magnetic field, 



cos 9 e sin k z 
vfi = v F {— — e ab H — -e c 
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The momentum distribution of zero-energy quasiparticles 
is simply given by 



Kpd(0, v f ) = Re 



1 



y/T+ P(v F ,iuj n 
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(20) 



III. RESULTS 



A. Zero-energy density of states(ZEDOS) 

First, we will discuss the density of states right at zero- 
energy. In Fig. 1, we show our results for the rj depen- 
dence of ZEDOS ^1^1 f or j = n a calculated with the 
BdG and the PD methods. ZEDOS under a field in the 
nodal direction(a = 0) is denoted as v n (Q) and in the 
anti- nodal direction(a = ^) as f a (0). Both v n (Q) and 
f a (0) are normalized to the normal-state value. 
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FIG. 1: The rj dependence of zero-energy density of states 
calculated with (a)the BdG method and (b)the PD method. 
The solid line shows ZEDOS for H || node(a = 0), i.e.i^l, 
while the dashed line shows ZEDOS for H || antinode(a = f ), 
i.e. i^M . ZEDOS is normalized to the normal-state value. 



Two results show a similar behavior that in a high 
magnetic field or when H is close to H c2 (lowry), v n (0) > 
f a (0), while in a low magnetic field (high 77), v a (§) > 
i/„(0). 

In order to understand this crossover behavior, let 
us see the momentum distribution of quasiparticles con- 
tributing to the zero-energy density of states. Figure 2 
shows the angle-resolved ZEDOS obtained with the BdG 
method for several values of rj. In a low field region (Fig. 
2(a)), it is found that narrow regions in Fermi surface are 
responsible for ZEDOS. In the case of H || node(a = 0), 
a sharp peak appears at 9 = 90° which corresponds to 
the nodal direction perpendicular to the magnetic field. 
In contrast, the nodal quasiparticles running parallel to 
the field (6 = 0°,180°) give much smaller contribution 
to ZEDOS. In the case of H j| antinodc(a = j), we can 
observe nodal peaks of the same height at 9 = 45° and 
135°. 

With increasing magnetic field (decreasing 77) (Fig. 
2(b)), nodal peaks become broader, and the contribu- 
tion from the core states increases by a large amount, 
commonly for both a. This is because the decrease of or- 
der parameter amplitude makes it easier for low-energy 
quasiparticles to propagate independent of a. However, 
as for the quasi-particles running parallel to the ficld(i.c. 
9 = 0°, 180°), we can see a big difference between the 
two cases. In the case of H || node, contribution from 
the field direction becomes larger with increasing field 
strength, while in the case of H || antinode, contribution 
from this direction remains small. 

When the field strength is increased (77 is decreased) 
furthur (Fig. 2(c)), nodal peaks become too broad to be 
identified. In most part of the Fermi surface, the angle- 
resolved ZEDOS recover the normal-state value. Never- 
theless, appreciable difference still exists for the quasi- 
particles running parallel to the field. 

We can understand the nature of this difference by us- 
ing the idea of quasi-classical trajectories. In the quasi- 
classical theory, a trajetory is determined from the Fermi 
velocity and an impact parameter. ZEDOS can be ob- 
tained by summing over contributions from those tra- 
jectories. Roughly saying, if the sign of order paramter 
changes on a trajectory, finite contribution to ZEDOS 
arises due to the formation of Andrccv bound states. 
Now let us consider the trajectories parallel to the field. 
In the case of H | node, quasiparticles propagating along 
such trajectories feel no superconducting gap due to the 
gap node irrespective of impact parameter. On the other 
hand, in the case of H || antinode, those quasiparticles 
feel finite and spatially uniform gap, because they prop- 
agate along the vortices and never cross vortex cores. 
(Only part of the quasiparticles run right through a vor- 
tex core and feel no superconducting gap. However, con- 
tributions from such trajectories are considerably small.) 
Hence, in the latter case, those quasiparticles are ham- 
pered by finite and uniform order parameter and cannot 
contribute to ZEDOS, no matter how small(high) order 
parameter amplitudc(magnetic field) is. This is why ZE- 
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FIG. 2: Angle-resolved ZEDOS (normalized to the normal- 
state value) calculated with the BdG method. 9 is measured 
from the field direction. (a)rj — 2.11, (b)r? = 1.05, and (c)?? = 
0.21 



DOS is suppressed in the case of H | antinodc in a high 
magnetic field region. 

Next, let us compare the above results with the angle- 
resolved ZEDOS derived from the PD method (Fig. 3) 
by integrating Eq.(20) with respect to k z . In most part 



of the Fermi surface, we observe much the same behav- 
ior as the BdG results. However, as for the quasiparti- 
cles running in the field direction(# = 0°, 180°), the PD 
method gives larger angle-resolved ZEDOS than that in 
BdG method. We attribute this difference to the pre- 
scription we made in solving the BdG equations, that 
is, a neglect of c-axis dispersion term. This prescription 
corresponds to considering only the trajectories parallel 
to the ab-plane, i.e. limiting the c-axis component of the 
Fermi velosity k z to 0. Hence, our prescription enhances 
the contribution of the quasiparticles propagating in the 
direction (6 = 0,k z = 0) compared with other contri- 
butions of k z ^ quasiparticles. As we have discussed 
above, it is those quasiparticles with k z = that sup- 
press ZEDOS in the case a = \. Therefore, we obtain 
smaller angle-resolved ZEDOS from the BdG method. 
For comparison, we plot the angle-resolved ZEDOS ob- 
tained from Eq. (20) with k z = in Fig. 4. Comparing 
it with Fig. 2(a), a quantitative agreement can be seen 
in 9 ~ 0° and 180°. 

Furthur, we comment on the difference in ZEDOS ob- 
served in a high magnetic field region^ <~ 0.1) between 
Figs. 1(a) and (b). Since the behavior of ZEDOS in this 
magnetic field region is dominated by the quasiparticles 
parallel to the field, our BdG analysis overestimates the 
difference in ZEDOS. 

Here we briefly summarize the main results in this 
subsection. In a low magnetic field region, ZEDOS is 
dominated by nodal quasiparticles which have a finite 
momentum normal to the field. In this region, we have 
f a (0) > v n {Q), since in the case of H || node, two of 
the four nodes are parallel to the field, thus unable to 
contribute to ZEDOS, while in the case of H || antin- 
ode, all the four nodes can contribute to ZEDOS. In a 
high magnetic field region, on the other hand, the differ- 
ence between ^ a (0) and v n {Q) comes from the behavior 
of quasiparticles running in the field direction. In the 
case of H || antinode, those quasiparticles are hampered 
by finite and uniform order parameter and cannot con- 
tribute to ZEDOS, while in the case of H || node, they 
can. Therefore, we have v n {Q) > ^ a (0) in this field region. 



B. Density of states(DOS) 

Next, we will show our results for the density of states 
at finite energy. DOS at finite energy is particularly im- 
portant when we discuss the experimental data, because 
it depends on the density of states at < e < ksT. In 
Figs. 5 and 6, we show DOS obtained with the BdG and 
the PD methods, respectively. 

Let us first discuss the high magnetic field region (Figs. 
5(b) and 6(b)), We can observe a sharp rise in DOS at 
|e| ~ Ao rcminicent of a coherence peak in the case of 
H || antinode, while not in the case of H || node. This 
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FIG. 3: Angle-resolved ZEDOS (normalized to the normal- 
state value) calculated with the PD method. 6 is measured 
from the field direction, (a)r; = 2.1 1 , (b)r/ = 1.05, and (cjr/ = 
0.21. 
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FIG. 4: Angle-resolved ZEDOS (normalized to the normal- 
state value) calculated with the PD method with k z = for 
7? = 2.11. 



character of DOS reflects the behavior of the quasipar- 
ticles running parallel to the field. Since they feel finite 
order parameter in the case of H |j antinode, they tend to 
form a coherence peak. This coherence-peak like struc- 
ture appears more clearly in the BdG result due to our 
prescription as discussed in the previous subsection. 

In this high magnetic field region, we observe v n (t) > 
u a (e) for < |e| < Ao independent of the calculational 
methods. Therefore, if the thermal conductivity or the 
specific heat is measured in this field region with a rotat- 
ing magnetic field, maximum will be observed when the 
field is applied in the nodal direction. 

The results for a low magnetic field region are shown 
in Fig. 5(a) and 6(a). In this field region, as we have 
discussed in the previous subsection, u a (0) > v n (0) at 
zero-energy. However, as shown in Fig. 5(a) and 6(a), 
u a (e) and ^„(e) cross at |e| ~ 0.2A . For |e| > 0.2A we 
have the opposite inequality v n {e) > v a {e). This means 
that when we take experimental data at a temparature 
T > 0.2Ao with a rotating magnetic field within the 
plane, we will observe very small angle variation, because 
the effects of lower-energy DOS and higher-energy DOS 
cancel each other. 



C. Interpretations of the experimental data 

In this subsection, we will discuss how to interpret the 
experimental data of the thermal conductivity and the 
specific heat, in connection with our analyses. 

So far the thermal conductivity^ and the specific 
healiiS of Sr2Ru04 have been measured under a rotat- 
ing in-plane field by several groups. In the experiments 
of the magnetothermal conductivity, they found no angle 
variation, except in the vicinity of H C 2. They attributed 
the angle variation near H C 2 to the anisotropy of H C 2 it- 
self, and denied the existence of vertical line nodes in 
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FIG. 5: Finite-energy DOS (normalized to the normal-state 
value) calculated with the BdG method. (a)»7 = 2.11 and 
(b)77 = 0.21. 



FIG. 6: Finite-energy DOS (normalized to the normal-state 
value) calculated with the PD method, (a)r; = 2.11 and 
(b)rj = 0.21. 



Sr2Ru04. On the other hand, in the experiment of the 
specific heat, fourfold oscillation is found at a lower field, 
in addition to the angle variation near H c2 . For the angle 
variation at the lower field, the specific heat shows max- 
imum at H| | [110] , while near H C 2 maximum is observed 
at H| | [100] . These measurements seem to provide incom- 
patible results. However, our theoretical analysis gives an 
answer to this contradiction in the following way. 

Here we discuss why the thermal conductivity does not 
show the fourfold oscillation in the low field region. As we 
show in the previous subsection, the reversal of v n (e) and 
u a (e) occurs in this field region at e ~ 0.2Ao- Thus, in 
order to observe the anisotropy of ZEDOS in this field re- 
gion, contribution from higher energy part of DOS must 
be removed. However, the thermal conductivity is mea- 
sured at rather high temperatures T = 0.2 — 0.3T C , where 
the effect of high energy part of DOS mixes inevitably as 
we note in the previous subsection. Hence, it is no wonder 
that the anisotropy of ZEDOS cannot be observed. On 
the other hand, the measurement of the specific heat was 
conducted at a lower temperature T < 0.1T C , where the 



specific heat is sensitive to ZEDOS. Furthurmore, note 
that the specific heat is more sensitive to the low energy 
part of DOS than the thermal conductivity. These will be 
the reason why they can observe the fourfold oscillation. 

Combining the above considerations, we can determine 
the position of line nodes. In the measurement of the 
specific heat, the maximum was found at H || [110], while 
from our calculations we have u a (0) > J-Vi(O) in a low 
magnetic field region. Therefore, we conclude that the 
line nodes exist in the direction [±100] and [0±10]. 

In sections IlIIA and HUB , we have shown that v n {e) > 
i^a(e) for < |e| < Ao in a high magnetic field region. 
Therefore, in this field region, there is a possibility that 
we can observe a fourfold oscillation both in the thermal 
conductivity and the specific heat, aside from the effect 
of the anisotropy in H C 2- 

In summary, we can explain the data of specific heat 
and thermal conductivity simultaneously by assuming 
the vertical lines nodes along [±100] axis. The line nodes 
in this direction is expected in the 7-band from the re- 
lation between the Fermi surface and crystal symmetry 
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(ReflU), and has been confirmed in microscopic analy- 
sis (Refjl^). Strictly speaking, these are not line nodes, 
because excitation gap is small but finite. However, tiny 
gaps serve as line nodes in the finite temperature. 

IV. CONCLUSIONS 

We studied the density of states in the vortex state of 
a layered superconductor with vertical line nodes on the 
cylindrical Fermi surface under a field parallel to the ab- 
plane. We investigated the angle variation of DOS with 
changing field strength. Bogoliubov-de Gennes equation 
and an approximate analytical method due to Pesch and 
Dahm were solved. We found that a field in the nodal 
direction gives larger zero-energy density of states in a 
higher magnetic field region, whereas a field in the anti- 
nodal direction results in larger zero-energy density of 
states in a lower magnetic field region. This crossover 
phenomenon is naturally understood in terms of the mo- 
mentum distribution of quasiparticlcs. In a higher field 
region, under a field applied in the anti-nodal direction, 
quasiparticles running parallel to the field is hampered 



by spatially uniform order parameter, and thus gives a 
smaller contribution to ZEDOS compared with the case 
in which the field is applied in other directions. In a 
lower field region, nodal quasiparticles not parallel to the 
field contribute to ZEDOS significantly. Therefore, when 
the field is applied in the anti-nodal direction, four such 
nodes are available, which leads to the larger ZEDOS 
compared with the case of H | node, where only two 
such nodes are available. 

We investigated the angle variation of the density of 
states at finite energy. We found that in a low magnetic 
field region, DOS at e > 0.2Ao shows maximum when the 
field is parallel to the node, while ZEDOS shows mini- 
mum for this field direction. On the basis of this fine 
structure of DOS, we discussed why the thermal conduc- 
tivity does not show fourfold oscillation in the low field 
region. Finally, combining the experimental data of the 
specific heat and our analyses, we conclude that Sr 2 Ru04 
has vertical line nodes in the direction of the a-axis and 
the b-axis. 

We are grateful to K. Deguchi for teaching us his ex- 
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discussions. 
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